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Abstract. 

The purpose of this contribution is to show that some of the basic ideas of turbulence can 
be addressed in a deterministic setting instead of introducing random realizations of the 
fluid. Weak limits of oscillating sequences of solutions are considered and along the same 
line the Wigner transform replaces the Kolmogorov definition of the spectra of turbulence. 
One of the main issue is to show that, at least in some cases, this weak limit is the solution 
of an equation with an extra diffusion (the name turbulent diffusion appears naturally) . In 
particular for a weak limit of solutions of the incompressible Euler equation (which is time 
reversible) such process would lead to the appearance of irreversibility. In the absence of 
proofs, following a program initiated by P. Lax [L], the diffusive property of the limit is 
analyzed, with the tools of Lax and Levermore [LL] or Jin Levermore and Mc Laughlin 
[JLM] , on the zero dispersion limit of the Korteweg-deVries equation and of the Non Linear 
Schrodinger equation. The three authors are extremely happy to have the opportunity to 
publish this contribution in a volume dedicated to Walter Strauss as a mark of friendship 
and admiration for his achievement. They hope that this paper concerned with non linear 
fluid mechanics, non linear instabilities and inverse scattering, will find its place in the 
different domains that have interested Walter. 
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1. Introduction. 



Two types of objects appear frequently in the theory of turbulence for a fluid defined in 
an open set O C R d , d = 1, 2 or 3 : Models of Turbulence and Turbulent spectra, the first 
ones is used in most of practical numerical simulations and the second proposed as a tool 
for the understanding of the phenomena. 

One of the most classical models of turbulence is the so called k — e model which reads : 

k 2 

d t U + V X (U®U) + V X P - isA x U - c u V x [—(V x U +* V X U)] = 0, 

V x -U = 



and 



d t k + UV x k - ^-\V X U + (V x Uy\ 2 - V x -[c u —V x k] + e = , 
d t e + UV x e - ~y\V x U + (V x Uf\ 2 - V x ■ [c 3 —V x e] + c 2 - = . 



(2) 



In the above system q denote several constants usually given by experimentation. Equation 

(1) is the standard incompressible Navier Stokes equation modified by the introduction of 
a "turbulent viscous" term : 

k 2 k 2 
uV x u + c u [—{V x U +* V X U)\ ; v T = v + c v — (3) 

where the space-time dependent scalar quantities k and e are defined by the coupled system 

(2) and interpreted in term of local energy fluctuation and local enstrophy. 
Intuitively the formula encompasses the idea that the fluid by its self interaction produces 
some averaging effect and therefore that the average U is a smooth quantity solution of an 
equation with a viscosity greater than the initial one. Therefore one conjectures and for 
some variant of the system (1) (2) (cf. [MP] and [LM]) proves the relation : 

v T > v . (4) 

This relation becomes important when v is very small. This situation corresponds to large 
Reynolds numbers. In some sense the idea is to have for the average an equation which 
would become valid when, due to the complexity, the initial equation cannot be computed. 
Such effect becomes crucial when v = 0. All what remains is the turbulent aspect of the 
viscosity. It has to be positive otherwise the turbulent model is an ill posed problem (in 
the sense of Hadamard) and the solution cannot be computed. Furthermore it gives one 
of the many examples where the limit of a reversible system (the incompressible Euler 
equation) becomes irreversible. 

There is up to now no mathematical analysis of the range of validity of such formulas. The 
" phenomeno logical" proof usually given introduces some randomness in the description of 
the motion of the fluid. The randomness plays an even more crucial role in the definition 
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of the spectra of turbulence as given below. A random family of realizations u(x,t,u>), is 
said to be homogeneous whenever the tensor 



(u(x + r,t) <g) u(x, t)) 

is independent of x and t and it is said to be isotropic whenever it depends only on \r\. 
The homogeneity hypothesis implies that the tensor 

R{k) = j ' u(x + r,t,u) <g> u{x,t,u)e~ ir ' k dudr (5) 

is independent of x and t. The isotropy hypothesis implies that it is given by the formula : 



R{k ' t) = 4^W {I -lkF ] - (6) 



The expressions (5) and (6) are starting points for most of the phenomenological theory 
of turbulence according to Kolmogorov [Ko] and Kraichnan [Kr] . 

2 Weak convergence and Wigner Measures. 

As said in the introduction one of the purpose of this contribution is to show that the 
recourse to randomness is not compulsory for the above construction. At variance, as 
already observed by Lax [L] , one could consider the weak limit of a deterministic sequence 
of oscillatory solutions of the fluid mechanics equations and use recent results concerning 
defect measures and Wigner transform (cf.[Ge], [LP] and [Ta]). 

The first observation is that with the homogeneity hypothesis, the formula (5) can also be 
written as : 



R(k,t) = / u(x + -,t,u) ® u(x - -,t,u)e %r ' k dudr 
J 2 2 



(7) 



and therefore the right hand side turns out to be the Wigner transform of u. As above it 
is a symmetric positive tensor. However the progress is that this tensor is can be written 
in term of a local quantity R(x, t, k, u) : 

R(x, t,k) = < R(x, t, k, lo) > 

^ f r t h 

R(x, t, k,u>) = / u(x H — , t, u>) <E> u(x ,t,u>)e lT ' dr 

J z z 

and homogeneity hypothesis can be relaxed. 

In fact this idea already appeared as a basic ingredient of the contribution of D.C. Besnard, 
F.H. Harlow, R.M. Rauenzahn and C. Zemach [BHRZ]. In spite of the fact that it is a very 
natural approach, to the best of our knowledge, it has not been used elsewhere in turbulence 
theory. Furthermore the Wigner transform, in space, or (after time localization) in time 
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is well defined for any solution u or sequence of solutions u n and one has the standard 
formulas : 

R n (x,t,k)= / u n (x+ -,t) ®u n {x- -,t,u)e~ ir ' k dr , 

(9) 



u n (x, t) <g> u n (x, t) = (-^) d / R n (.i\tJv)(lk 



R e n (x, t,r) = J e(s)u n (x, t + |) <g> u n (x, t - |, oo)e- lST ds, 9 e P(R+), 6»(0) = 1 , 
u n (x,t) ®u n (x,t) = (-^) y R e n (x,t,r)dT . 

The two last formulas of (9) are important because they involve instead of space corre- 
lations, time correlations which are the quantities more commonly involved in practical 
experiments. 

The only physical a priori estimate, uniform with respect to the Reynolds number, is the 
energy estimate. 

\ I \u n (x,t)\ 2 dx < \ I \u n (0,t)\ 2 dx < C. (10) 
Jn Jo. 

This observation is valid in particular for the 3d Euler or Navier Stokes equation as proven 
by Di Perna and Lions (cf [Li] section 4.3). It is also valid even in 2d for the Navier Stokes 
equation with viscosity going to zero, when the natural viscous boundary condition u = 
on <9Q is assumed. In particular Grenier has constructed solutions to the Navier Stokes 
equation in the half plane O = R y+ x R x with vorticity blowing up in the L°° norm when 
the viscosity goes to zero([Gre] theorem 2.1.). As observed by Grenier this is an instability 
phenomenon for the Prandtl layer which is of the same type as the non linear instability 
of the Euler equation proven by Friedlander Strauss and Vishik [FS V] . 
With the estimate (10) one concludes that, up to the extraction of a subsequence, u n 
converges in weak*L 00 (H t+ , L 2 (0)) to a limit u. However due to the above considerations, 
in many cases, one will have : 

lim u n (x, t) <E> u n (x, t) = u(x, t) ® u(x, t) + R tur b(x, t) (11) 

n^oo 

with R tur b(x,t) 7^ 0. In a follow up of ideas of Peter Lax [L] one could consider that the 
appearance of the tensor Rturb plays the role of the Reynolds stress tensor as the emergence 
of turbulence in a deterministic approach. In fact it is the Defect Measure of the sequence 
u n . The tensor Rturb is symmetric positive definite and one has : 

Rturb{x, t) = lim (u n -u)<S> (u n - u) , (12) 

n^oo 

or, with the Wigner transform of u n — u, 

Rturb(x,t,k) = lim / (u n - u)(x + ^,t) <g> (u n - u)(x - T -,t)e~ %r ' k dk (13) 
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and 

lim u n (x, t) (g) u n (x, t) = u(x, t) (g) u(x, t) + ( — ) d / Rt U rb(x,t,k)dk . (14) 

The formula (14) displays the natural link, for our purpose, between the defect measure 
and the Wigner transform. Observe that the left hand side of (13) is the natural local and 
deterministic "avatar" of the Kolmogorv spectra for random turbulence and it is natural 
to conjecture that it will inherit the basic properties of isotropy and scale law for the 
dependence in k. The isotropy hypothesis is made plausible 3d by the following remarks : 
A necessary condition for R tur b(x,t) to be non zero is that curl u n becomes unbounded 
in the neighborhood of (x, t) then : 

(i) for non zero viscosity it has been shown by Constantine and Fefferman [CF] that it is 
much more the oscillations in direction of the vorticity than its size that are responsible 
for instabilities in the fluid, 

(ii) exterior constant Coriolis force stabilizes the fluid when the Rossby number goes to 
infinity as shown by Babin Nicolaenko and Mahalov [BNM] and others. Therefore a de- 
composition of the vorticity according to the formula : 

curl u n = O n + Q n (15) 

with O n having a constant direction and a modulus going to oo while u n remains bounded 
should not be possible. 

In the sequel of this section we consider in the two dimensional case, sequences of solutions 
to the Euler equation with the impermeability condition : 

d t u n + V x {u n <g> u n ) + V x p n = , 

(lb) 

n ■ u = on 90, V x • u n = , 

or of the Navier Stokes equation in a domain O of R 2 with the viscous boundary condition : 

d t u n + V x (u n <g> u n ) + X7 x p n - V n A x ll n = 0, ^ 
«„ = on 90, V x • u n = . 

In both cases existence and uniqueness of such solutions are well established facts with the 
hypothesis 

u n (x,0) = u° n (x) e L 2 (0), W x -u° n = 0,u;° = V x < e L°°(0) , (18) 
in the first case and with the assumption 

u n (x, 0) = u (x) e L 2 (0), V x -w = , 
in the second case {y n > 0). 

However in these two cases the turbulent Reynolds tensor may be present in the limit : 
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(i) if one considers a sequence u n of solutions of the Euler equation with initial data 
uniformly bounded in L 2 (0) but with initial vorticity unbounded in L°°, 

(ii) in the second case if, as already mentioned above, one keeps the initial data fixed but 
let the viscosity v n goes to zero. 

In both cases the limit satisfies the equations 

d t u + V x (u ®u)+ V x R tU rb + V x p = 0, V x -u = , (19) 

with 

Rturb(x,t) = lim (u n (x,t) - u(x,t)) <S> (u n (x,t) - u(x,t)) . (20) 

n— >oo 

Introducing the trace : 

T 

(19) is changed into 



pll I p22 

^turb ^turb 



d t u + X7 x (u®u) + V x S tU rb + V X P = 0, V x -u = . (21) 

with Sturb denoting a tracefree tensor and P = p + T. The space of tracefree tensors is 
of dimension 2 and assuming that that the limit u is a smooth function, this space has a 
natural basis given by the matrix : 

J(V« + V*«)=f lm 9 ^ , Ud* 2 ui + d Xl u 2 )\ 
2V ' \2{dx 2 ui + d Xl U2) d X2 u 2 J 

and an orthogonal complement 

Hu)= (W^ + d Xl u 2 ) d X2 ui \ 

V d X2 u 2 -^{d X2 u 1 + d Xl u 2 ) J 

Therefore there exist two space-time depending functions vturb( x it) and S(x,t) such that 
one has : 

Sturb = Vturb (VU + V*u) + 8${u) (24) 

and the equation (21) becomes the equation : 

d t u + V x (u®u) + V x (u turb (Vu + V*u)) +V x (5$(w)) + V x P = 0, V x -u = 0. (25) 

As discussed above, the mechanism of "creation of turbulence" should be isotropic and 
this would imply that the tensor Sturb is invariant under Galilean transformations and 
therefore proportional to (Vw + V*it) reducing (25) to a diffusive type equation : 

d t u + V x (u®u) + V x (u turh (Vu + V*u)) +V X P = 0, V x -u = 0. (26) 

Observe however that in any case, assuming that on the boundary u is zero whenever the 
vturb is positive one has the energy estimate : 

\d t \ \u(x, t\ 2 dx + / vturb{ x i t) | V x u\ 2 dx = . (27) 
Jq Jo. 



A turbulent model would be obtained by coupling equations (25) with a system of equations 
which would determine the function vturb- The necessary condition to obtain a such well 
posed system is that the function vturb(x,t) is non negative, a property which does not 
follow for the fact that the tensor Rturb is itself non negative. 

In the case of the Euler equation, with a sequence of initial data having unbounded 
vorticity one would derive an irreversible problem as the limit of reversible equations. 
Eventually from the formula (27) one deduces the 



Proposition 1 In the above configurations, assume that the sequence of initial data 
Un(x) converges strongly to uq(x) in L 2 (Q), then for T > 0, the following assertions are 
equivalent : 

(i) the sequence u n {x,i) converges strongly in L 2 (0 x [0, T]), 

(ii) vturb{x->t) is identically zero on O x [0, T], 

(iii) one has : 

T 

I / v tU rb(x,t)\V x u\ 2 dxdt <0. (28) 
Jo Jn 

Proof The only non classical point is the fact that (iii) implies the strong convergence or 
equivalent ly that one has : 

lim inf i / / \u n (x,t)\ 2 dxdt = \ [ [ \u(x,t)\ 2 dxdt . (29) 
n ^°° Jo Jn Jo Jn 

With the classical energy estimate (for any given n) and the relation (27) one deduces the 
inequalities : 



— / \uo(x)\ 2 dx > lim inf h [ [ \u n (x,t)\ 2 dxdt = h [ [ \u(x,t)\ 2 dxdt . 
2 Jn n ^°° Jo Jti Jo Jn 

> 7T / \u (x)\ 2 dx- [ [ i>turb{x,t)\V x u\ 2 dxdt 
Jvl Jo Jq 



(30) 



and (29) follows from (28). 



3 Positivity versus non positivity of the Diffusion Coefficient for the small 
Dispersion limit of KDV and NLS flows. 



In the absence of a systematic theory it seems worth while to study, following the program 
of Lax, the issue of the positivity of the turbulent coefficient on the dispersive limit of the 
KDV and NLS equations using explicit formulas given by the inverse scattering theory. 
For the KDV flow one considers the problem : 

u\ — 6u c u e x + e 2 u e xxx = 0, with initial data u e (x, 0) = uq(x) (31) 
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and the for the NLS flow the problem : 



e 2 

ieu e t + -u e xx + (1 - \u e \ 2 )u e = 0, 

Six) 

with initial data u e {x, 0) = A(x)exp(i ) . 



(32) 



With the introduction of the functions : 

12 - — ~ ie 



«T - 1, and p e = —{u e u% - u x u e ) , (33) 



the NLS equation is equivalent to the system : 

p\ + p\ =0, 



ft + i^ + C) =^(P^ogp% x ) x , 



(34) 



The equations (31) and (32) are time reversible and (34) is a reversible perturbation of the 
the usual isentropic compressible Euler equation. For e going to zero the functions u e , p e 
and p e converge weakly and the following notations are introduced : 

u = weak lim u e , 
u 2 = weak lim (u e ) 2 , 



p = weak lim p e , 
a = weak lim u e , 
u e2 o e2 

Q(p^p1 = (^- + V)' 

p z 



(35) 



and one obtains the equation : 



Q = weak lim Q(p e , p e ) 



U 2 s ,-u 2 -u 2 



- ea x ( y ) - ea x ( y - -) = o (36) 



and the system 

Pt + fix =0 , 

— 2 —2 

<9t£ + <9 X (^- + ^r) + d x (Q- Q(p, p)) =0 . 
p 2 

In the region where strong convergence occurs one has 



(37) 



or 

Q-Q(p,fl)) = 0. (39) 

As expected in these regions (36) is (up to a simple change in the x variable ) the Burgers 
equation and (37) is the compressible Euler equation for isentropic fluids. 
On the other hand it is known (cf [LL] and [JLM]) that the strong convergence does 
not hold everywhere. The regions where strong convergence fails are called the Whitham 
regions. By a simple convexity argument, (observe that both the functions 

u — > u 2 and (p, p) — > Q(p, p) 

are convex), one has in the Whitham region : 

2 -2 

y - y > and Q - Q(p,p)) > 0. (40) 

To analyze the possibility of the appearance of "turbulent viscosity." one writes (36) and 
(37) in the following form : 

y2 (u) 2 

u t - 3((w) 2 ) x - d x (u turb d x u) = 0, with v twb (x,t) = — ^— — (41) 

and 

Pt + fix = , 

p 2 p 2 

Pt + {—T + ~^)x — QxiyturbPx) = , 

Q-Q(p,ji)) 

With Vturb\X, t) = — . 

d x p 

The existence of any kind of turbulent model requires that v tU rb be non negative which is 
equivalent here to the property : 

d x u{x, t)>0, (43) 
in the Witham region for the KDV dispersive limit and : 

Q-p^))> 0> (44) 
d x p 

in the Witham region for the dispersive NLS limit. This is also equivalent to 

d x p > . (45) 

Since diffusion properties may appear on a larger time scale it is natural to explore the 
properties (43) and (45) for large time. 

Such program is done below using the tools of the inverse scattering following [LL] and the 
conclusion will be the fact that such properties are satisfied depends on the initial data. 
The starting point are the following theorems : 
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Theorem 2. ([LL]) Let u{x,t;e) solve 

u t - 6uu x + e 2 u xxx = 0, (46) 

with initial data u(x,0;e) = uo(x) belonging to the Schwartz class, strictly negative, and 
with only one minimum point xo, at which Uo(xo) = — 1. Let x±(rj) be defined for < rj < 1 
by 

uo(x±(r])) = —rj 2 and X- < xq < x + . (47) 

Define the function 

r x+(ri) 



m = L M RiRP fc ' (48) 



for < rj < 1. Then, 
(i) the weak limit 



u(x,t) = lim u(x, t; e) 



exists. 

(ii) As t goes to infinity, for x such 5 < x/t < 4 — 5, with 5 is any given small positive 
constant one has 

u{x,t) = -^{{^ 2 ) + o{l/t). (49) 

(iii) As t goes to infinity, fora;/t<0orx/t>4 one has 

u = 0(t~ 2 ) . (50) 

This theorem is stated and proved in the third paper of Lax and Levermore [LL] pages 
810-815. Furthermore, with some conjecture on the uniform effect of the remote part of the 
initial data on the solution of KDV equation the authors adapt their asymptotic analysis 
to the initial data : 

/ x f — 1 if x < 0; /ri \ 

M0(X) = \0 iix>0. (51) 

which correspond to the shock profile for the Burgers equation. They obtain for the weak 
limit the formula : 

u(x, t) = — 1 for x < —6t , 

u(x, t) =s{-) for - 6t < x < At , (52) 
u(x, t) = for At < x , 

where the function £ i— > s(£) can be computed in term of complete elliptic integrals. 

Explicit numerical computations done in [LL] on the formula for s(£) indicate that this is 

an increasing function on the interval [—6,4]. 

The theorem 2 has a counterpart for the NLS dispersive limit. 

Theorem 3. Let u(x, t; e) solve the NLS flow : 

ieu t (x, t; e) + — -u xx (x, t; e) + (1 - t; e) \ 2 )u(x, t; e) = , (53) 
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Six) 

u e (x,0) = \u Q {x)\ 2 exp(i^-) , (54) 

with 

p (x) = l + \u (x)\ 2 , (55) 

and 

p (x) = d x S(x) , (56) 

belonging to the Schwartz class. Let also assume, for simplicity, that the initial data are 
"single well" in the following sense (cf [JLM]). Introduce the functions 

r±(x) = \d x S±A{x), (57) 

and assume that r_ has only one maximum X max while r+ has only one minimum \ m in 
with the relation : 

-1 < r-(x) < X max < Xmi n < r + (x) < 1 . (58) 
Define the numbers x ± (X) according to the formula : 

for - 1 < A < A max ,r_(x±(A)) = A, x_ < x + , 
for \ min < A < 1, r+(x±(X)) = A, X- < x + . 

Then, the weak limit 

p[x, t) = lim e ^o\u(x, t; e) | 2 + 1 
exists. Furthermore in the Whitham region 

x 

— £ ( — 1, A max ) U (Xmirn 1) , (59) 



as t — > oo, one has : 



4 

p(x,t) = 1 <f>(x/t)(l - (x/t) 2 ) 1 ' 2 , where 

Tit 

, m _ r +W A-l/2(r + ( g ) + r_( g )) 
^ J "A.(A) (A-r + ( S ))^(A-r_( S ))V2 ds ' 



(60) 



and elsewhere one has, p ~ 1. 

Proof : The proof is given in [K] and for sake of completeness the main steps are recalled 
here. The existence of the weak limit is proved in [JLM] . Following a method suggested in 
[LL] one begins with the multisolitons formula for fixed e and then let e go to zero. 
For fixed e, the long time behavior of \u\ 2 is as follows [FT, pp. 168-176]. In the solitonless 
region \x/t\ > 1 or X max < x/t < A m i n , one has 

\ u \ 2 = i + o(r 1 ' 2 ). (6i) 
li 



As t — > oo. In the Whitham region, the solution is a multisolitons solution : 

\u(x,t;e)\ 2 ~ 1 - E^{s(x - r] n t - x n ,rj n ), where 

- I-?? 2 (62) 

COS^CCl-T/ 2 ) 1 ^^) ' 

with exponentially small error. The rj^s are the associated eigenvalues (of the underlying 
Dirac operator) and the x n 's are some phase constants of no importance. 
The width of each soliton is 

°(jT^fW y (63) 

By Weyl's law for the distribution of eigenvalues in (—1, X max ) U (A m j n , 1) as e — > 0, 

- Vn = —7-T , (64) 

nvn) 

where f/ n G (r] n ,r] n+1 ). 

Peaks of solitons are located at r\ n t. As t — > oo, they are separated by 



net 

<P(Vn) 



(65) 



so for large t they are well separated. 

The wave number 77 of the soliton that peaks at x at time t is 77 = x/t, if t is large and 
either -1 < x/t < A ma:E or X min < x/t < 1. 
Therefore the density of the solitons is 

» (66) 
net V ; 

The area between a soliton and the line it = 1 is 

4e(l - r/ 2 ) 1 / 2 ~ 4e(l - (x/t) 2 ) 1/2 , (67) 
so the asymptotic area density which is given by 1 — p is the product of (66) and (67) : 

1^/^(1 -(x/t) 2 ) 1 / 2 . (68) 

7T 

Hence, the asymptotic formula for the weak limit p follows. 



^From the above statement several observations can be made concerning the appearance 
of a positive turbulent viscosity in the limit equation satisfied by u for the KDV dispersive 
limit and by (p, fi) for the NLS dispersive limit. Such positivity would be related to the 
appearance of irreversibility in a weak limit of reversible models. 
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As said above the construction of [LL] (section 7 page 817) shows that a shock profile as 
initial data produces in the limit a smooth solution with a "turbulent viscosity". 
Following Theorem 2, one considers initial data u(x, 0; e) = uq(x) belonging to the Schwartz 
class, strictly negative, and with only one minimum point xo, at which Uq(xq) = — 1. For 
large t the solution is asymptotic to 

-T-<t>((^) l/2 ) with <j>(rj)= / — " — ^-r^dx. (69) 

Ant^Ht' ' ^" J x _ {v) (-uo(x)- 77 2 )V2 ^ > 

Therefore one has : 

d x u~ L_0'(^_). (70) 

yr(4t)2 (4t)2 

And the "turbulent diffusion" hypothesis requires that 

(j>'( V ) <0,V V e [0,1]. (71) 

With uq(x) = — e~l x l' 3 one has : 

/•( 21 <)^ „ 

^ »/ = 2 / 7 — ? r 2 U /2 & • ( 72 ) 

On the table 1 the values of 4>{rf) and for the following (3 exponents : (3 = 1, 3/2, 2, and 4 
and with 9 steps r\ = /clO -1 , 1 < k < 9 

In the first case is decreasing, in the second case its variation changes, and then is 
increasing in the two last cases. 

For the NLS flow one considers initial data satisfying the hypothesis of Theorem 3 and 
observes that the existence for large time of a diffusive regime would be given by d x \x > 
in the Witham region and with the conservation law : 

Pt + fix = , 

this is equivalent to pt < in the same region. With Theorem 3 this condition is equivalent, 
for large t to the relation 

<9 A (A0(A)(1-A 2 )^) >0,VAe [-1,+1]. (73) 

with 0(A) given by 

m ~L m (A-r +W )'/ 2 (A- r -W)'/ 2 rfs ' (7 "> 

Solution with initial data having zero momentum and a symmetric density with an frac- 
tional exponential rate of convergence at infinity are analyzed : 

p G (x) = l-\e-\ x \\ S(x)=0. (75) 
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The fact that the momentum is zero gives r + (A) + ?"-(A) = 0. The x symmetry of the 
density remains true for all time and all e therefore it is enough to consider the behavior 
of 

(A0(A)(l-A 2 )i)0(A),A=|, 

for 

< \ ml n < A = | < 1 . (76) 

One has : 

! p (-log 2(1- A)) 13 A 2 ,/Tl - X 2 ) 

A0(A)(l-A 2 )2)0(A) = 2 y o X* - (l - l e -\x\y dX ■ (77) 

The numerical computation given on the table 2 are done for (3 = 1.5, 2 ,3 and 3.5 with 
A = x/t varying from \ m i n = 0.5 to 0.9 with step 0.1. They indicate that the "turbulent 
regime" appears for (3 = 1.5, (3 = 2, (3 = 3 but does not hold for (3 = 3.5. 



(3 = 1 


= 3/2 


(3 = 2 


(3 = 4 


5.88251 


2.19345 


1.75226 


.78565 


5.47775 


3.65183 


2.07024 


1.09422 


5.06441 


3.37626 


2.29727 


1.36578 


4.63711 


3.09141 


2.47717 


1.62689 


4.18879 


2.79252 


2.62703 


1.89375 


3.70918 


2.47278 


2.75579 


2.18372 


3.18159 


2.12106 


2.86860 


2.52387 


2.57400 


1.71600 


2.96898 


2.97357 


1.80404 


1.20273 


3.05934 


3.73515 



Table 1 Numerical computation for the dispersive KDV limit. 

Values of 4>(rf) are computed for (3 = 1, (3 = 3/2, which appears as a critical case, (3 = 2 
and (3 = 4 with rj = lO^k, 1 < k < 9 



(3 = 1.5 


1.10957 
1.62180 
2.16222 
2.69501 



(3 = 2 


1.29635 
1.64627 
1.97843 
2.21615 



(3 = 3 


1.48400 
1.63429 
1.76762 
1.76986 



(3 = 3.5 
.0 

1.53471 
1.62240 
1.70143 
1.64733 



Table 2 Numerical computation for the dispersive NLS limit. 

Values of A^>(A) are computed for (3 = 1.5 (3 = 2, (3 = 3 and (3 = 3.5 which appears as a 
critical case with A = 0.5, 0.6, 0.7, 0.8, 0.9. 



Conclusion. 
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In this contribution it has been shown that some of the basic questions of the statistical 
theory of turbulence could be formulated in a deterministic setting with the introduction 
of sequence of weakly converging solutions. The counterpart of the turbulent spectra being 
the Wigner transform and the turbulent diffusion being related to defect measures. Explicit 
computations done on integrable classical integrable system indicate that for these models 
it is not always possible to construct a "formal" turbulent equation. At this point of our 
analysis it depends on the behavior of the initial data and in particular on the fact that 
they should not be too much concentrated (their decay for |x| going to infinity has to be 
not too small) . It is worth while to notice that with the convenient conjectures of [LL] the 
shock profile leads always to a diffusive regime. 
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